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ABSTRACT 


It is shown how to use a reliability shorthand to get 
analytical results for small systems. The underlying proba- 
bility distribution is the exponential function. An algorithm 
1s described to convolve exponentially distributed random 
variables, and two computer programs are given to obtain 
numerical results for convolutions and mixed convolutions. 
After describing systems in shorthand notation, the programs 
can be used to compute reliabilities for those systems. 
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RENNE ODUET TON 


Every product eventually fails and, sooner or later, is 
no longer usable for its original purpose. The rate at which 
this deterioration process proceeds is called the failure rate. 
One is interested in predicting the time to failure for 8 
component or an entire system. It is impossible to predict 
the time to failure precisely, but one can model this process, 
viewing the time to failure as a random variable which has 
a probability distribution that can be specified. There 
exist numerous life distributions which provide realistic 
probabilistic descriptions of the behavior of components or 
whole systems. 

This paper is concerned only with the exponential distri- 
bution which plays a fundamental role in reliability. This 
distribution has many properties which facilitate the modelling 
process, and so encourage its application in many fields. 
Where not directly applicable, it is often used as a bound. 

According to the exponential distribution, a component 
has a constant failure rate. An equivalent property is that 
the remaining life of the component is independent of its age: 
the "memoryless" or "no aging" or "no wear" property. This 
means that a used component is always as good as a new one. 

Convolving life distributions for components to get the 


final life distribution for the system is rather difficult, 





even for the exponential distribution. For an increasing 
number of components it's no longer possible to compute a 
closed form life distribution. This paper will show how to 
use a reliability shorthand notation to get analytical results 
for small basic systems, and later use an algorithm to com- 
pute numerical results with an Apple II+ microcomputer. The 


Pascal programs are provided in the appendices. 





II. RELIABILITY SHORTHAND 


A. NOTATION AND DEFINITIONS 
The shorthand notation "EXP [A1] t EXP [A4] " denotes the 
distribution for a random variable Tj+T,. T, has an exponen- 


tial distribution with failure rate እ ገ and T, has an exponen- 


mer distribution with failure rate iA... Ty and T 


2 are assumed 


2 
to be independent. 
In the following T always denotes a random life length 


and 


Probability (Tot), = 1 - F(t)‏ - اشام 


is the survival function as a function of time. 
The convolution of life distributions corresponds to 
the summation of components with independent random lives 


to get the life distribution of a system. 


10 





EM E 1 شك‎ 1  ፣ FOR A ONE COMPONENT SYSTEM WITH SPARES 


A. GENERAL DESCRIPTION 

There are two ways for a simple system, consisting of 
one active component and one spare, to survive a specified 
time t. Either the active component survives time t, in 
which case the spare need not be used, or the active com- 
rent tabls at some time 5 (0 < s < t) and the spare takes 
over and doesn't fail in the interval from s to t and thus 
completes the mission. 


ነ FÐ) and F_(t) be the exponential survival functions 


2 
for component one and two respectively, and f, (t) and f£ (t) 
their respective densities. Then the survival function F(t) 


for the system is denoted by: 


(t-s)f,(s)ds 


E 
H= P (t) + E 55 


where the integrand shows the probability that component one 
fails at time s (£, (s)) and component two takes over for 

at least the length of time t-s. The following shows some 
examples how to get the survival function in closed form for 


amibterent basic convolutions. 


L] 





Boe BASIC SURVIVAL FUNCTIONS 


LS sten ben Neidentical Components 


The shorthand notation for this system is: 


EXPO + EXP (AJ) +>... + EXP [À] 


Whenever an active component fails a spare takes over, as 
long as there is a spare. When the n-th component fails, the 


whole system stops functioning. 


lor n — 2: 
መ — t 
F(t) = F (t) d 4 ጀጋ (55) ፻1 (5) 65 
-ìt € OTA (t-s), As 
= e + f e le “ds 
0 
= (l+atje**, Cac 


Mier density function for this survival function is: 


d © = D E atte" 
= =[(1+እቲ) (-እፀ ^7) + SS 
EE E 
"Bs 
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Kor N = 3: 
The survival function for the two component system 
will be convolved with the survival function of the next 
component to obtain the survival function for the three 


component system. 





E 
Hom rs "ሪራ 
0 
which leads to: 
PIS x: 
F(t) = (l+ÀXt+—" e 7 ! O 


Bor all n: 
Continuing in this manner one gets the general solution 


Bor the survival function of this kind of system: 


n 

) (At) 
" = 
mie) = O ` t 

Bee 


i-l 


| V 
O 


This class of convolutions always gives the well-known 


Belang[n,A] survival function. 


2. System with N Dissimilar Components 


The shorthand notation for this system is denoted 


by: 


EXP [11] T EXP [45] +... EXP [A] 
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For n = 2: 


The survival function for the active component is 


The survival function for this system is 


| V 


I V 


E 
Hr MEI cT Ets). (syds 
0 
x t ipee (t-s) =A,S 
Zu. Pa. he Tas 
1. 
0 
st À -\_t È “(እ>ለ.)8 
cc ou AS f (A.-A.)e EE ds 
À. “እ 2 
1 2 
=A,t A = A E (AL-AL)t 
e ee Å L 0 2 se Em S 
1 2 
zo t -A.t -入 .七 
1 2 1 
1 Ay እጋ)ወ + 1,e i.e ) 
(Ay = እጋ) 
>ለ ጋኒ =A,t 
(Aie is እጋ 5 
= E O 
(5 — A.) M = 
1 2 
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For n = 3: 
ik Ar-suleltġ Survival function for the two component 
system will be convolved with the survival function of the 
next component to obtain the survival function for the three 
component system. 


The result is 


B AN -A.t N -入 .tt 
Í m uc EORR kd (À er 233) © I 
sc n 2፡1 ST 
. Ay A Ag 
AA) A TA) 


Bør all n: 
Continuing in this manner one gets the general solution 


for the survival function of this kind of system: 


፲ 


PA me 
gc P 1 
j#ži 


The two general survival functions for systems built from 
components which all have dissimilar failure rates or all 
have the same failure rate are very useful for computations 
Eastonctaonsot systems which are formed from sub- 
systems of the above two kinds. 

It will be shown in the following section how to use 
the above results to get a closed form for the survival 
function of a system which can not be computed with only one 


of those general formulas. 
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3. Systems With Identical and Dissimilar Components 


Shorthand notation: 
EXP [41] - EXP [41] 十 EXP [A] ፦ EXP [1] 


Let the random variables T be distributed in the following 
way: 


Ti - EXP [11], T, - EXP [A1], 


T4 = EXP [AJ], T4 E EXP [A]. 


There are different approaches available for computing the 
Survival function for this system. 

One can compute the survival function for the first 
two and then the second two components, using the general 
formula of Section III.B.l and then convolve those two sur- 
vival functions; or one can compute the survival function of 
the first two components and then add the other components 
one at a time. The latter approach will be used here to 
give a general idea of the procedure, and to lead to an 
algorithm for the numerical evaluation of the reliability of 
systems consisting of so many components that a closed form 
is no longer attainable or not attainable in an efficient 
manner. 

The first two components (፲1 = Ti, T1 = T t T) 


2 1 


have as shown in III.B.1 the survival function 


16 





+ Ë 
| + | Er. (t-s)fr,(s)ds 
i 55 ر‎ 
-àt 
(811 ተ-51 3) ; to ao (1) 
wach: 
511 = 1 and 5813 = እገ 
Adding the next component leads to T; = Ti + T, + T, 
ama the survival function 
F 3 Lm 
የ = ም 
T3{t) = Ep, (t) E e c s) £4, (s)ds 
2 0 3 2 
-入 -七 CE Ns) -ÀS 
te + + fe 2 12se Í as 
1 0 il 
which leads after the integration to 
2 2 2 2 
= ፦ - 一 人 
E (A, 22 A2I FAJA) እገኣጋ)ር Àt Ax P 
Emp + (€) = > = Me 
3 (እ1”ላ3ጋ) (^i እጋ) 
or in the form of equation (1) to 
_ =A E “335 
Bp ም = 811151) 十 521 
with 
22 2 
ኝ _ እ 72),A, ۴ 2142-2 _ እ ገ 
ll ` D 12 7 Á 2 = _ 27 
(417A 2) (Að 9) (3.1 Àg) 
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Adding the last component leads to T, = Tj+T,+T,+T, and 


the survival function: 


Ë 
E. (t) pM uU (t=s) f, (s) ds 
T4 = o Ta = 


A,t “እጋኒ 
(a, ta, tie 十 5319 


E “A, (t-s) -À.S 


1 
bi y S Speci o eS 


“ላ 38 
- እ 35318 ] ds 


Integrating and collecting terms gives the final survival 


function for the whole system: 


E “እገኒ “ላጋኒ 
E = (a, ¡ta¡,t)e + (831532) P t > 0 
with 
2 2 
l i 入 - 5 l À jÀ 9 
ME pa E * 12 ` 2,37 
(A) = à 3) (እ2”* 1) 
3 2 2 
_ እገ ” 2እ1ላ3 ሺ i ATA, 
I oa NM 2 
(44734) (Ay እጋ) 


This section demonstrates that one can get the survival 
function by adding one random variable at a time and thus 
computing the new survival function. Each summation then 


gives a new or updated ZEE Even though the integration is 
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not shown explicitly, one can imagine that this method of 
convolving survival functions is very tedious and it is very 
easy to make mistakes. So the next section introduces an 
algorithm to compute all the aij” which then makes it easy 


to evaluate the survival function at any time. 


C. ALGORITHM TO COMPUTE THE CONVOLUTION OF EXPONENTIAL LIVES 


The shorthand notation to be solved is 


EXP [1] 4 EXP [1 I + ,.. ተ EXP [A1] (ni terms) 
+ 5፳፻ [3.3] + EXP [1 „ ] + ,.. ተ EXP [3.3] (n, terms) 
DNT ا‎ ie 
+ EXP [À L] + EXP [2 ] HIE. EXP [à] (nj terms) 


The general solution for the survival function is of the form 


i -A,t IE ME 
F(t) = A, (t)e + A, (t)e To... + A, (t)e puto O 
with 
2 3 m 
A (t) = a] + ajot + agt + 514: + ,.. ተ jin, f ; 
A, (t) = 831 十 533 一 53 3” 十 534 + ,.. ተ ns ; 
-1 
A, (t) = a + a t + a ው +a t3 +... +a CU E 
K ki k2 k3 k4 kn, 


An example shows how to get the survival function for the 


following shorthand notation: 
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EXP [Ay ] + EXP [41] 十 EXP [Àj] 
+ EXP [X] + 5፳፻ [3.3] 


十 EXP [À] 


Survival Function: 


mit) = (a tta..t je 22003 


113112: 


Ger 0 。 


An algorithm to compute all the factors 511/ jor ---። «31 
and thus make it possible to convolve all the different 
exponential functions was earlier introduced by Gursel [Ref. 
1], and is given here in a slightly modified form. This 
algorithm adds one exponential random variable in each run. 
The above example can be solved in this way by looping six 
times through this algorithm. 

Netation for the algorithm: 

K = number of dissimilar failure rates (3 in above 


example); 


l. = i-th type failure rate, 1 = 1,...,Kk; 


nn. = number of random variables having the 9-th 
type failure rate; 


Lx = failure rate of the ie-th lambda, ie = l,...,K 
n, = current number of random variables having the 
i-th failure rate, n, = 1,...,nni; 
a = j-th coefficient of the i-th polynomial 
NEN Minde for current polynomial and currently 


entering life. 
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Algorithm to compute the coefficients: 


Initlalization: 


ሠ 
ዘ 


O 2552 TK where jJ = 12... ያምፅ 


Jik os MA 
l 
n, = 0, for all 1 where i = l,...,K 
TID ut : እ; for all i where i = l,...,K 
nn, forall where 1 = l,...,K 
After the first run: 511 > Ty n, = 1 
eae, Until n. = nn. 
ie ie ie ie 
1. Update the coefficients a. EOL KE susen. 
le,k ie 
= * = == መዉ 
e irs aan 1), where m Nie” for 
j = 0,... n, 3 
2. Update the coefficients a. 
le ,1 
n n ` 
_ E aiid; i E 
ارد د د كل‎ ` Pl j ` 
1 4 izie i ie j=2 5) 
i=1 ifie 


3. Update the other coefficients a.. for all i,k where 
K ፕ | x O 


for 1 = | ere ar 
Ru. uH 


1 
Sia. = S ee for all i # ie, 
£ I 
. = a. À. -ma. À. =À, for all 1 where 
irn, | mul r9 im R E i Z ie and n. l 


for m = nj-]， 
ገ 一 1,...,(n,-=1) 
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D. PROGRAM FOR NUMERICAL EVALUATION OF THE RELIABILITY 
FOR A ONE COMPONENT SYSTEM WITH SPARES 


The above algorithm is used to compute the reliability of 
a one component system with similar and/or dissimilar spares 
for any given time. The input to the program, which is 
written in PASCAL and listed in Appendix A, will be explained 
here with the same example as above: 

EXP [A1] Sr EXP LA, I i EXP [1 I 


b EXP[A,] 3 EXP [A5] 


+ EXP[A3] 
Let's assume that እገ = 0.1, እጋ 10727 Az = 0.3 and we want 
دن‎ |1169 the reliability for t = 0, 10, 20, ..., 90 hours. The 


prompts of the program are typed in large letters. 


"HOW MANY DISSIMILAR LAMBDAS DO YOU WANT TO USE?" 

3 

"WRITE DOWN ALL THE DISSIMILAR LAMBDAS WITH AT LEAST ONE 
BLANK IN BETWEEN (INPUT NUMBERS < 1 WITH A LEADING ZERO) :" 
1 0.2 0.3 

"WRITE DOWN THE NUMBER OF EACH OF THE LAMBDAS WITH AT LEAST 
ONE BLANK IN BETWEEN:" 

1 د 

"HOW MANY TIME POINTS DO YOU WANT TO BE COMPUTED? (<=20) :" 
10 

"AT WHAT TIMES DO YOU WANT THE SURVIVAL FUNCTION TO BE 
EVALUATED? (TIME SHOULD HAVE THE SAME DIMENSION AS THE 


LAMBDAS AND SHOULD BE AT MOST 20 INTEGERS BETWEEN 0 AND 70. 
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WEE UPPER LIMIT 70 IS REQUIRED ONLY, IF A GRAPHICAL 
PRESENTATION IS DESIRED LATER) :" 

0 10 20 30 40 50 60 70 80 90 

"DO YOU WANT A SIMULATION TO BE DONE PARALLEL TO THE 
TOMNEUTATIONS? (SERIOUS SIMULATION NEEDS A LOT OF TIME!) HIT 
or "YES" AND "0" FOR NO:* 

l 

"HOW MANY NUMBERS DO YOU WANT TO USE AS BASIS FOR THE 
SIMULATION? :" 


1000 


After the computations the program will produce the following 


ere pu : 


FOLLOWING LAMBDAS WERE ENTERED: 0.1 0.2 0.3 


+ OF EACH LAMBDA: 3 2 b 

Results: 

RELIABILITY ( T» 0) - 1.000000 (SIM = 1.000) 
BEBTABILITV(T > 10 ) = 0.995799 (SIM = 0.996) 
EMBTABILITV(T 20) = 0.925259 (SIM = 0.925) 
RELIABILITY ( T > 30 ) = 0.739626 (SIM = 0.728) 
RELIABILITY ( T > 40 ) = 0.508910 (SIM = 0.489) 
DENS ApBrLITY (T » 50 ) = 0.311271 (SIM = 0.334) 
RELIABILITY ( T > 60 ) = 0.174439 (SIM = 0.173) 
BEETABILITV(T > 70 ) =. 0.091597 (SIM = 0.068) 
RELIABILITY ( T > 80 ) = 0.045784 (SIM = 0.036) 
RELIABILITY ( T > 90 ) = 0.022028 (SIM = 0.020) 


"DO YOU ALSO WANT A GRPAHICAL PRESENTATION OF THE DATA? (HIT 


"m 1 " FOR "YES " AND " 0 " FOR "NO te ) e n 
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Mae graphical output will not be shown here, but will be 


generated from the program if desired. 


l. 


Some notes on the program: 


The program is set up to allow 20 dissimilar lambdas. 
It is necessary to type leading zeros in PASCAL if one 
inputs decimal numbers. 

The simulation generates n numbers (in this example 

n = 1000) for every random variable, to produce n 
observations on the whole system. 


Because of 


the behavior of the exponential function, 

- the sequence of the lambdas in the input, 

- the method of adding always only one random variable 
at a time, 

- the possible closeness of the values of the lambdas 
and 

- the possibly large number of the dissimilar lambdas, 

the program might have to compute with very large and 

small numbers and could produce erroneous results. 

Therefore a simulation is provided within the program 

to see if the computations are close to the simulation. 

Fixed limits which avoid errors can't be given, because 

all the possible errors are closely related. Fewer 


dissimilar lambdas for example would allow that their 
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values could be closer together, or that less care is 
needed with the sequence in the input, etc. 
But after a long test phase and a comparison with results, 
produced with the IBM 3033, some hints can be given: 

1. The input sequence of the lambdas should go from the 
smallest to the largest, 

2. There should be no more than 20 lambdas altogether, 

3. The difference in value of the dissimilar lambdas should 


be larger than 0.05. 
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0 00 RELIABILITY FOR A SYSTEM OF ACTIVE COMPONENTS AND SPARES 


A. MIXTURE OF CONVOLUTIONS 

Before solving systems with active components and spares, 
the mixture of convolutions will be explained briefly. Let's 
assume that we want to compute the reliability of a parallel 


System with two different components, like that in Figure 14. 





Figure 1. Simple Parallel System 


The system starts with a life distributed as EXPÍA tro]. Then 


1 


there are two different possibilities: either component fl 


fails with probability P. = 和 1/ (入 1]+ 入 7) and the residual life 


1 
is distributed as EXP[\.]. This is the upper branch of 


Figure 2. Or component #2 fails with probability 


P- = እጋ/ (A tA.) and the residual life is distributed as 


2 
EXP [۸,1۰ 
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EXP [A] 


EXP [1 j +A. J - 





EXP [A] 


Figure 2. Shorthand Notation for Parallel System 


Using the "MIX"-notation, the life can be written in the 


shorthand notation as: 
* * 
EXP [A,*314] * MIXÍP,*EXP[A,] , P *EXP [41] ). 


Because of the distributive law, which is valid for these 


mixtures, one can transform this equation to 
* * 十 
MIX{P,*EXP[A.] + EXP[A,+A,]),P,*(EXP[A,] + EXP[\, X,1)) 


The first and the second expressions in parentheses multiplied 
bv their respective probabilities, is a mixture of convolu- 
tions. The following section demonstrates how those mixtures 
can be solved numericallv even if an analvtical solution is 


not feasible. 
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ES DESCRIPTION OF THE SYSTEM 

This section is concerned with computing the reliability 
of a system of components. The different components might 
be supported by spares, which have the same or a different 
failure rate as the original component. The following example 
will be used to illustrate the idea of describing the system 
in the reliability shorthand notation. The listing of a 
PASCAL program for these types of systems is provided in 
Appendix B. The program incorporates the method of program 
#1 to add the components, one at a time, to get first the 
convolution of the exponential distributions and then via 
the mixture the reliability of the whole system. 

Example: 


The system consists of two parallel components (#1 and 


#2) in series with a third component (#3). A spare is pro- 
vided for component +1. The components have the failure rate 
入 A, and A, respectively. The spare has the same failure 


I 2 3 


rate as the component it can replace. 





Figure 3. Combined Serial/Parallel System with Spare 
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C. SHORTHAND BRANCHING DIAGRAM 


À 
1 
EM | (1) 
እ1ተላለጋትላጋ“ DEE 
^2 
— _EXP[\_+_] (2) 
E 2a 3 
s EXP [A +A +A] + : 
— B 
À .+À +À —— ZERO (3) 
pag 入 ] + 入 2+A3 
= 
Å ነ i Ss (4) 
2 
——— T 
EXP [Aj +A, +À 31+ US ES [A ተእጋ]ተ 
J Oo 3 Å» 
一 一 一 一 ZERO (5) 
0 ከ. 
A FA tA 
— 
Figure 4. Shorthand Branching Diagram for Serial/ 


Parallel System with Spare | 


The numbers in the numbered parentheses in the figure 
above describe the different endpoints (leaves) in this nota- 
tion. The connections between the original system (root) 
and the leaves will be called branches. In the following: 
the connection from the root to leaf #1 is called branch #1, 
from the root to leaf #2 branch #2, etc. 

The original system has a life distributed as EXP [A,+A,+A31] . 
Three different things could happen now. Either component 


tl fails with probability Ay ZA] FA FA 3) or component #2 fails 
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Pech Probability እጋ/ (3:1 ተላ ጋተጳ 3) or component #3 fails with 

probability Az/ (ARAS +A 2) Let's continue in branch #1 and 
assume component #1 failed with the above probability. Now 
the spare can replace the failed component and the system's 


life is again distributed as EXP [A] +À +^3] . The system is now 


2 
in the same state as at the beginning and the same arguments 
hold. Let's assume new component fl fails with probability 


እገ/ AG +A tA), then the remaining life is distributed as 


2 
EXP [AHA]. If component #2 fails (the probability of this 
failure is hy lA tantra) ), the remaining life is distributed 
as EXP [Aj +À q]. The third possibility, component #3 fails with 


probability እጋ/ NA FÅ th), brings the failure of the whole 


2 
system, denoted by the "ZERO"-notation. This says that there 
is a remaining life of zero. The reason why this third branch 
isn't simply omitted is that even if there is no remaining 
life, there is a specified probability that this event occurs. 
And this has to be taken into account to compute the total 
life of the system. 


The other three branches of the above figure are developed 


analogously. 


D. THE SURVIVAL FUNCTION 
Let's denote the survival function of branch #1 by F, (t), 
and the probability that the events in branch #1 occur in 


this special sequence as P Then the life of branch #1 is 


1 


distributed as P *(EXP[)A +A +À,] + EXP [à] +À 


1 +A 2] + EXPIA,TA,1) 


2 


30 





BMOMNE survival function of this system is the following 


mixture: 
= = m km m km = = 
F (t) MIX{p, Fj (t),p, OS F4(t),p, Fj (€) PG FF (€) , pc *Fc (t) ] 


For the computation of the coefficients for every single 
survival function the algorithm of Section III.C can be used 


again. 


E. PROGRAM FOR THE NUMERICAL EVALUATION OF THE RELIABILITY 
FOR SYSTEMS WITH SEVERAL ACTIVE COMPONENTS AND SPARES 


This section describes the program for computing the 
reliability for systems as described above. While reading 
this part, one should look at Figure 4. The program is 


interactive. Only the input of Figure 4 will be explained 


here. 
Prompt Input 
"How many branches?" 6 CR (Carriage return) 


EC Branch : 


l. Element CEE SCR O CR 

Ze Element : IMERMITERM2TER 3 CR O- CR 
3. Element : Re FERIER VER 

4. Element : Mere RT CR 3 CR OCR 
5. Element : لت ل‎ OR CR O CR 

6. Element : UCR 

e Branch : 

1, Element : CRIS CR O CR 

2. Element : EEE? CR 37 CR O CR 
3. Element : Mer 2 2 312 0 CR 
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Element 
Element : 


Element 


Branch 2 

Element : 
Element : 
Element : 
Element : 
Element : 


Element : 


Branch 
Element : 
Element 
Element 
Element 
Element 


Element : 


Branch 
Element 
Element 
Element 
Element : 
Element 


Element : 


Branch : 
Element : 
Element : 
Element : 


Element 


> ፦፦ ፦ Wo 
No FB IM ثم‎ 
W N W N 


O FP Ñ P F Hp 
N HF ሠ ፦ ሠ 
O O 


NTE Co FE 


O PON FP 0 O 


O H w o 





The branches in the figure count from top to bottom and the 
elements in each branch from the leaf to the root. The first 


element in the first branch is EXP [À „+À } and one inputs only 


3 
the two indices of the lambdas. The second element in the 


same branch is Ay IAS FÅ +h.) and one inputs first the index 


2 
of the lambda in the numerator and then the indices of the 
denominator. After each element input a "0". If there is no 
sixth element as in branch fl to branch #5, no fourth element 
as in branch #6, or no element as the ZERO elements in 
branches #3, #5, and #6, Simply input a "0". Because of the 
data structure used in the PASCAL program, the system always 
has to be split up into pieces, so that there is only a single 
failure in each step. Because of the algorithm used, there 


፡ከ በቢ ገ 'ነ፦ be more than 20 elements in each branch to avoid 


ENDOIS. 
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V. SUMMARY 


The reliability shorthand is a convenient way to describe 
the life of different systems. This shorthand then is the 
underlying data structure to compute numerically the relia- 
BEY for systems for which it is difficult or impossible 


to get the survival function in closed form. 
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APPENDIX A 


PROGRAM CONVOLUTION; 
كان‎ TRANSCEND, APPLESTUFF; 
CONST 
ASTERISK=” 1” 
BLANK 
DASH 
CINE 
TYRE 
ARR 
ARI 
VAR 
MAEREKRED ARRAYE1..20,1..20123 OF REAL; 
NRO. REL, RELSIM :ARR; 
NINIT,NI,T:ARI; 
TO IE, K, TMAX., SIMMAX, Y, X. COUNT. VAL: INTEGER; 
03007 1 J.EACT.NN.NKK,NNI.CHGR,CHSIM: INTEGER; 
ED :INTEGER; (X INDEX FOR T-ARRAV ¥) 
SØM RELI:FEAL: 
PRINTVAR: CHAR; 
(¥ L = ARRAY OF DISSIMILAR FAILURE RATES x) 
(% NICTI = NUMBER OF LIVES HAVING THE ITH FAILURE RATE 9 
GCI. J] = COEFFICIENT FOR THE POLYNOMIAL x) 
(ሄ T = ARRAY OF TIMES AT WHICH THE RELIABILITY WILL BE 
BBNMPUTED 8) 


? 


ae. A9 MB “Mes 


> 
2 ? 
a 


1 2 


PACKED ARRAYE1..203 OF REAL; 
PACKED ARRAYL1..20] OF INTEGER: 


PROCEDURE ATTENTION; 
VAR 
DURATION, 1: INTEGER; 
MESAS ARRAYE1..103 OF INTEGER; 


BEGIN 
ENT CHIL1J:- 15; 
ELCJ:= 17; 


IMTIEHESJ:- 19; 
BIPFEHLE4J:72 20; 
BĦAMENEOSJ:—- 22; 
DURATION:= 30; 
FOR I:=1 TD 5 DO 
NOTE(CPITCHLI3,DURATION) ; 
END; 


PROCEDURE DECIDEZ: 

BEGIN 
WRITELN:WRITELN: 
WRITELN(7'DO VOU WANT TO CONTINUE WITH THE SAME PROGRAM? 
re "لا"‎ 2’); 
WRITELN(”?DO YOU WANT TO COMPUTE MIXED CONVOLUT IONS? 
TRES Ve Gre SE 
ETP 1): 
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END; 


PROCEDURE INPUT; 
BEGIN 
WRITELN; WRITELN; 
IDEROX Y (0,12); 
WRITELN(”HOW MANY DISSIMILAR LAMBDAS DO YOU WANT TO USE? "5; 
READLN (K) ; 
WRITELN; 
WRITELN(' WRITE DOWN ALL THE DISSIMILAR LAMBDAS WITH AT 
LEAST ONE BLANK IN BETWEEN: ’); 
WRITELN(” (NUMBERS < 1 NEED A LEADING ZERO!)?’);3 
EUR 13-1 TO K DO 
READ(LEIJ2); 
WRITELN; 
WRITE C WRITE DOWN THE NUMBER OF EACH OF THE LAMBDAS WITH AT 
LEAST ONE BLANK IN’); 
WRITELN(’ BETWEEN: 7) 5 
FOR I:-1 TO K DO 
READ(NINITLII): 


WRITELN; 
WRITE( HOW MANV TIME POINTS DO YOU WANT TO BE COMPUTED? 
"ጠይ 20) 5 ``); 


READLN(TMAX) ; 
WRITELN(”AT WHAT TIME DO YOU WANT THE SURVIVAL FUNCTION TO 
BE EVALUATED?”); 
MMMIETEEN(T (INTEGER TIMES BETWEEN 0-50 IN THE SAME UNITS AS 
THE LAMBDAS) ?); 
FOR I:=1 TO TMAX DO 
BERD(TEI)); 
WRITELN; 
WRITELN(’DO YOU WANT A SIMULATION TO BE DONE PARALLEL TO 
THE COMPUTATIONS?"); 
WRITELN(” (SGERIDUS SIMULATION NEEDS A LOT OF TIME')')3 
do: FOR YES OR "o" FOR NO : au 
READLN(CHSIM) s 
WRITELN; 
MEWEHSIM 7. 1 THEN 
BEGIN 
WRITELN (7? HOW MANY NUMBERS DO YOU WANT TO USE AS BASIS 
FOR THE SIMULATION?”); 
READLN (SIMMAX) ; 
END: 
WRITELN: WRITELN; 
END: 


ERBEEBDUÜRE OUTPUT; 
BEGIN 
AE Fr OLE OWING '".k:5,'" LAMBDAS WERE ENTERED :?”); 
FOR Y:=1 TO K DO 
BETVETEEI]:5:2); 
WRITELN; 
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# OF EACH LAMBDA :"); 


WRITE (’ 
FOR I:=1 TO K DO 


WRITE(NINITLIJ:5); 
WRITELN; WRITELN; WRITELN: 


end; 


PROCEDURE INITIALIZE; 


BEGIN 
KKs=13 
WHILE kk<=k DO 

BEGIN 
ALKK,13:=1.0; 
NICKKIJ:=1;5 
JJ:-1; 

WHILE JJ<=K DO 
BEGIN 
IF JJ“ >KK THEN 
BEGIN 
ከ .1 ነ52=ሣ1.‹.., JF .ህህ3/ (.ኒቪፒህህጋ3ገ3-ኒፒኮር]3)2 
END: 
JJ:=JJ+1; 
END 
KE:=ERF I; 
END; 
Kke=kK; 
IE:=1; 
END: 
PROCEDURE SUMMATION; 
BEGIN 

SUM: 20.0; 

memels=1 TO F DO 
BEGIN 

IT 1> IE THEN 
BEGIN 


ግ. SUMSALI.IJECLIJ/(QLLIJ-ELIEJ3)5; 
IF NILI3?-2 THEN 


BEGIN 
FACT:=1; 
NKKs:-NILIJ; 
FOR Il:=2 TO NKK DO 


BEGIN 
FACT:-FACT£(11-1); 


RER ELIE) THEN 
5DUMSSSUMTAPI,IIJXEACTEXLEIEJ/EXPCIIXLN( 


CrO- CCTTE J?) 
El SE 


BEGIN 
IF (IT MOD 2)=0 THEN 
SUM: =SUM+ALI , II JXFACTXLCIEJ/EXP ( 


JIXLN(ABS(L[1I-L[1E3))) 
BESE 
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SUM ባባ ዘገ+፡)(1. 1118 ዕር: ቆ3ጊ! 1..1/ (-ክር.፦ሥ( 
IIXLN(ABS(LLIJ-LLIEJ)))) 


END; 
END: 
END; 
END; 
END: 
pate. i ji:=SUM+AL IE, 13: 
END; 


ዚዜ ከ CEDURE SUR1; 
BEGIN 
“ስሙ 
WHILE I«zk DO 
BEGIN 
IF I<>1E THEN 
BEGIN 


A AL ITONICIASSECCIED/ (LETES—-LOEII); 


TEEN ]IEI)JZ>1 THEN 
BEGIN 
J:=1; 
NNI:-NILIl-J: 
WHILE J<=(NICTII-1) DO 
BEGIN 


ALI.NNIJ:-(LDIEJXALI,NNIJ-NNIEAEI,NNI-11)/ 


EN An EN 
3. -=J)+1: 
NNI:=NILI]-J; 
END; 
END; 
END: 
is-Tel: 
END; 
END; 


PROCEDURE ONEATATIME; 
BEGIN 
IF NINITLIEJZNILIEJ THEN 
IE:=IE+1 
ELSE 
BEGIN 
DI CTET: =NILC[IE]J]+1:; 
J: =O; 
REPEAT 
NNs=NICIEI-J; 
ALIE.NN3:-zLELIEJXALIE,NN-13/(NN-1); 
s=J+1; 
UNTIL NN=2; 
SUMMAT ION; 
SURI? 
END; 
END; 
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PROCEDURE RAND; 

VAR 
I, II, III, COUNT, J: INTEGER; 
R, TSIM:REAL; 


BEGIN 
፲ ፲ ፲ 2=ርን2 
COUNT : 20; 
REPEAT 
TSIM: =0.0;3 
CD 
FOR J:=1 TO K DO 
BEGIN 
FOR I:=1 TO NINITLJI DO 
BEGIN 
R:=RANDOM/MAXINT; 
ው ን ከ ፡=!ር 1 ግተ(>ሙ(1ህ‹( 2)/1.:01173))2 
END; 
113-1115 
END; 


Pee 11] 1111 1 THEN 
COUNT : =COUNT +1; 
11ل‎ 3-21 1 1+1 : 
UNTIL 1 1 1-5 1: 
BEESSIMLUITJ:-COUNT/SIMMAX; 
END: 


K UCEDURE SUB2; 
BEGIN 
BEEI:=0.0; 
IF TLITI70 THEN 
RELI:21.0 
ELSE 
BEGIN 
FOR I:=1 TO K DO 
BEGIN 
SUM: 20.0; 
NKK: ZNIL[ IJ; 
FOR J:=1 TO NKK DO 
SUM: =SUM+AL 1, JIXEXP ( (I-1)ALN(TC1ITI);S 
Pree 2s —=SUMZEXFP(-LEITIZTCITI); 
END; 
FOR I:=1 TO K DO 
RELI:=RELI+PROLI IJ; 
END; 
ከ - اع‎ ቃ1.ቄህ THEN 
RELCITIJ:=1.0 
ELSE 
REL ET: =RELT: 
IF CHSIM=1 THEN 
RAND; 
END; 
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PROCEDURE YAXIS; 
BEGIN 
IF (X=0) THEN PRINTVAR:-zL INE; 
IF ((X20) AND (Y=20)) AND (TL1J=0) THEN PRINTVAR: =ASTERISK; 
IF (X=-3) AND (Y MOD 2 = 0) THEN PRINTVAR: =’ 0? 3 
IF (X=-3) AND (Y=20) THEN PRINTVAR:=7’ 17’; 
IF (X=-2) AND (Y MOD 2=0) THEN PRINTVAR:=’.73 
IF (X=-1) AND ((Y=20) OR (Y=0)) THEN PRINTVAR: =’0O’ 5 
IF (X=-1) AND (Y=18) THEN PRINTVAR:="9"; 
IF (X=-1) AND (Y=16) THEN PRINTVAR: =’ 8’ 
IF (X=-1) AND (Y=14) THEN PRINTVAR:=" 7? 
IF (X=-1) AND (Y=12) THEN PRINTVAR: =’ 6’ 
IF (X=-1) AND (Y=10) THEN PRINTVAR:="5" 
IF (X=-1) AND (Y=8) THEN PRINTVAR:="4"; 
IF (X=-1) AND (Y=6) THEN PRINTVAR:="2"; 
IF (X=-1) AND (Y=4) THEN PRINTVAR:="2" 
IF (X=-1) AND (Y=2) THEN PRINTVAR:=”" 1”; 
END: 


‘ae 
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PROCEDURE GRAPHICS; 


REGIN 
COUNT :=13 
FOR Y:=20 DOWNTO 1 DO 
PEGIN 
SUR X:=-5 TD 76 DO 
BEGIN 
IF COUNT<=TMAX THEN 
BEGIN 
VAL: =ROUND (2O*RELE COUNT I): 
IF (VAL=Y) AND (X=TCCOUNTI) THEN 
BEGIN 
PRINTVAR: =ASTERISK: 
COUNT: =COUNT +1; 
END 
EL SE 
PRINTVAR:=BLANF; 
END 
ELSE 
PRINTVAR:=BLANK; 
YAXIS; 
WRITE (PRINTVAR) : 
END; 
END; 
Y:=03 
FOR X:=-3 TO 76 DO 
BEGIN 


IF (X MOD S =0) THEN PRINTVAR:=LINE; 
IF COUNT<=TMAX THEN 
BEGIN 
VAL : ZROUND (2OXREL LCOUNT 31) ; 
IF (VAL=Y) AND (X=TCCOUNTI) THEN 
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BEGIN 
PRINTVAR: =ASTERISK; 
COUNT : COUNT *1; 
END 
ELSE 
IF (X MOD 5 <> 0) THEN 
PRINTVAR: =DASH; 
END 
ELSE 
IF (X MOD S <> 0) THEN 
PRINTVAR: =DASHs 
IF (X=-3) THEN PRINTVAR:=?” 0” 
IF (X=-2) THEN PRINTVAR:=”.? 
IF (X=-1) THEN PRINTVAR: =’ 0’ 
WRITE (PRINTVAR) ፤ 
END; 
WPRITECC ር) 5 10 15 (0 e SL 二 二 40 45 
50 55 Sem: 
WRITELN(' Mee"): 
END: 


ae 40 a8 


BEGIN 
INEUT; 
OUTPUT; 
FOR IT:=1 TO TMAX DO 
BEGIN 
INITIALIZE; 
REPEAT; 
ONEATATIME; 
UNTIL IE=K+1; 
SIS 
+1 1. 1101 THEN 
mine RELIABILITY ( T >’,TEITI:3,* D 
Á S= ምርር! >1 ን ፤፣1 122254 ) ? 
ELSE 
RITE RELIABILITY TCT 2” TLITIL,? ) 


MEET BE 


> RELEIT2:8:4); 


END; 
WRITELN:WRITELN; 
ATTENTION; 
WRITELN(’DO YOU ALSO WANT A GRAPHICAL PRESENTATION OF 
BUE DATA’) >: 
CAREC HIT "1" FOR YES OR "0" FOR NO : E 
READLN(CHGR) 3 
WRITELN; 
WRITELN; 
IF CHGR=1 THEN 
GRAPHICS: 
DECIDE2; 
END. 
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APPENDIX B 


PROGRAM MIXED CONVOLUTION; 
USES TRANSCEND, APPLESTUFF 


CONST 
ASTERISK=? 4? ; 
BLANK =” ?; 
DASH = = le 
LINE =? !? 35 
MAX =10; 
TYPE 
ARR = PACKED ARRAYIT1..MAX] OF REAL: 
ARI = PACKED ARRAYS 1..MAX] OF INTEGER; 


DARR = PACKED ARRAYL1..MAX.1..MAXJ OF REAL; 

VAR 
AR:PACKED ARRAYIG. .MAX,0O..MAX,0O..MAX] OF INTEGER; 
à: PACKED ARRAYL1..MaàX.1..MAXJ) OF REAL; 
SUMREL:PACKED ARRAYLO..MAX3 OF REAL; 
N:PACKED ARRAYLO..MAX1 OF INTEGER; 
NINIT,N1,T, ZA, INP:ARI; 
PROB, RELTOT, LAMB: DARR:; 
MURO. REL, P2ARR>5 
BR. BRANCHES., COUNT, KK, IL,KKK, JJ, IEs INTEGER: 
11 . 1. FA, NN, NEE, NNI, CHGR, TMAX, Y, X. VAL: INTEGER; 
SUM. RELI, LAM, DEN: REAL; 
PRINTVAR, CH: CHAR; 
(€ L = ARRAY OF ALL LAMBDA *) 
(X PROLIJ = PROBABILITY OF EACH BRANCH 1) 
(¥ T = ARRAY OF ALL TIMES AT WHICH THE RELIABILITY WILL 
Se COMPUTED ¥) 
(x BRANCHES = NUMBER OF BRANCHES IN SYSTEM 4) 


PROCEDURE INPUT; 


BEGIN 
WRITELN(* HOW MANY DISSIMILAR LAMBDAS ARE IN THE SYSTEM? °’); 
READLN (K) ; 


WRITELN(? WRITE DOWN ALL THE LAMBDAS WITH AT LEAST ONE BLANK 
IN BETWEEN: ’); 
WRITELN(? (NUMBERS < 1 NEED A LEADING ZERO!)7); 
FOR I:=1 TO K DO 

READ(LEII); 
WRITELN(' HOW MANY TIME POINTS DO YOU WANT TO BE COMPUTED? 
( «2 20 ) : °); 
READLN(TMAX) ; 
WRITELN('AT WHAT TIME DO YOU WANT THE SURVIVAL FUNCTION TO 
BE EVALUATED?’ ); 
WRITELN(" (INTEGER TIMES BETWEEN ዕ-ሬዕ IN THE SAME UNITS AS 
THE LAMBDAS)7); 
FOR I:=1 TO TMAX DO 

READ(T[ 13): 
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WRITELN; WRITELN; 
END; 


PROCEDURE ATTENTION; 
VAR 
DURATION, 1: INTEGER; 
PITCH: ARRAYL1..103 OF INTEGER; 


BEGIN 
PITCHL13:-2 15; 
ESSICHE21:-7 17; 
EESPIHES1I:— 19; 


መ 11 ተ ባ12= 20: 
E CHCS]:= 22; 
DURATION: = 30; 
ከ [15 1!:2=1 TO 5 DU 
NOTE (PITCHLC II, DURATION) ; 
END; 


PROCEDURE OUTPUT1:; 
BEGIN 
WRITE(' FOLLOWING ',K:5,7 LAMBDAS WERE ENTERED 
Perm 1:-1 TO K DO 
MIDE CLIO: E): 
WRITELN; 
END; 


ERBBEDURE OUTPUT2: 
BEGIN 
FOR I:=1 TO BRANCHES DO 
BEGIN 
FOR J:=1 TO MAX DO 
BEGIN 
FOR K:=1 TO MAX DO 
BEGIN 
IF (ARLI„J.K2=0) AND (K7>1) THEN 
ARIT.J,K]l:==1; 
END; 
END; 
END; 
END; 


PROCEDURE GETDATA; 
BEGIN 
FOR I:=0 TO MAX DO 
FOR J:=0 TO MAX DO 
FOR K:=0 TO MAX DO 
ART1.J.KJ:=-1; 
NLOJ1:-71; 
WRITELN; 
WRITELN( 7 HOW MANY BRANCHES?” ) ; 
READLN (BRANCHES) 3; 
:=0; 
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REPEAT 
I:2I-*1; 
WRITELN(I:3,7. BRANCH: °); 
: =O: 
REPEAT 
J:=J+1; 
PI IELN(J:Z, . ELEMENT: 7); 
:=0; 
REPEAT 
K:=K+1; 
READ CZALK1) 3 
ARLI.J,K1:2Z2AEEK1; 
NLJI:TFK3 
UNTIL ARII.J.K 3203; 
UNTIL ((ARCI,J-I,NLJ-11170) AND (ARLI,J.NLEJJI120)) 
AND (ZAL11=0)3; 
UNTIL I=BRANCHES; 
QUTPUT2 
END; 


PROCEDURE COMPPROB; 
BEGIN 
FOR I:=1 TO MAX DO 
FOR J:=1 TO MAX DO 
PROBLI,JI:z-1; 


BEGIN 
DEN: =O; 
Kı=2; 
COUNT: =1; 
WHILE ARSI,J,KJ>O DO 
BEGIN 
WHILE ARLI,J,KIJO DO 
BEGIN 
DEN: =DEN+LCARCI,J,KI14; 
Ks=kK+15 
END; 
PROBCI.COUNTI:=LCARCI,J,113/DEN; 
COUNT : =COUNT +1 5 
J: =J+2; 


K: =2; 


END; 
FOR I:=1 TO RRANCHES DO 
BEGIN 
፡=1: 


PCIJ:=1; 
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WHILE PROBLI,J3>0 DO 

BEGIN 
PEIJTPITITPROBEI,;J3; 
32=ኒህተ12 

END ; 

END; 
WRITELN; 
END; 


PROCEDURE COMPLAMB; 
BEGIN 
FOR l:=1 TO MAX DO 
FOR J:=1 TO MAX DO 
LAMBLI,JJ:=-1; 
Mat) 5 
:=1; 
WHILE I<=BRANCHES DO 
BEGIN 
LAM:=0; 
K:-21; 
COUNT:=1; 
WHILE ARLI,J,K2>=0 DO 
BEGIN 
WHILE ARCI,J,KIS0 DO 
EEGIN 
LAM:=LAM+LIARLI,J.KJ23; 
r: =F+1; 

END; 
LAMBI I, COUNT I: =LAM; 
COUNT : =COUNT+15 
J:=J+2; 
k:=1; 

LAM: =O; 
END; 
15-1 15 
2=12 
END; 
END; 


ERBEEDURE SORT; 
VAR 
EMP: REAL; 
SORTED: BOOLEAN; 
BEGIN 
REPEAT 
و‎ 
SORTED:=TRUE; 
WHILE LLI+13>0 DO 
BEGIN 
ርን. ]› [1[ተ11 THEN 
BEGIN 
በ፲! :2=! [17): 
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:10 2:13 ]|2 : 13 1 !ا 
LEI+11: =TEMP;‏ 
SORTED: =FALSE;‏ 
END;‏ 
SIC;‏ 
END:‏ 
UNTIL SORTED=TRUE;‏ 
mee ት! ]1:=፳-1:‏ 
mut21:-2-1;‏ 
END:‏ 


PROCEDURE INITIALIZE; 
BEGIN 
KK:=1; 
WHILE KK<=k DO 
BEGIN 
ALKK,11:71.0; 
NI[KkKI:=1; 
JJ: =1; 
WHILE JJ<=k DO 
BEGIN 
IF JIS >KK THEN 
BEGIN 
ከ .፡፡፣?:>ሙን፣ኮ ርጋ 11311 ህህ1ሪ (.!ህህ- [ርኢፌ1)25 
END; 
JJ:=JJ+1: 
ENDS 
pr SKr FI: 
END; 
KK: =l; 
15 5 -1 : 
END; 


PROCEDURE SUMMATION; 
BEGIN 
SUM: =O. O; 
FOR is=1 TO ር DO 
BEGIN 
ገ IZZIE THEN 
BEGIN 
SUM: =SUM+ACL 1, 1718#[(:111/ (11 ] ]-'[15ጅ]3): 
IF NICIIJ>=2 THEN 
BEGIN 
FA:=1; 
NKK: =NILII: 
FOR IIs72 TO NKK DO 
BEGIN 
FAs=FAR(II-1)3 
IE LIIJSLLIE] THEN 
SUM: SUMĦALI,IIIXFAXLLIEJ/EXP(TIŻLN( 
Lid deified) > 
ELSE 
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BEGIN 

IF (II MOD 2)=0 THEN 
SUM: =SUM+AL I „ I I JXFAXLLIEJ/EXP ( 
IIXLN(ABS(LLIJ-LLIEJ))) 

ELSE 
SUM:ZSUM*TALI,IIJEXFA*LLIEJ/(-EXP( 
IIXLN(ABS(LLIJ-LLIEJ)))) 

END; 
END; 
END; 
END: 
END: 
AE): =SUM+AL TE, 1); 
END; 


REDES EDURE SUBI: 
BEGIN 
1:21; 
WHILE I<=K DO 
BEGIN 
BE IIE TREN 
BEGIN 
BEILNRMEIJJ:-ORI,NILIJJELLIEJ/((LE IE2-LTL I2); 
IF NICTI< >i THEN 
BEGIN 
ل‎ : 21 : 
NNIS=NICTI—-J; 
WHILE J<=(NIL12-1) DO 
BEGIN 
ALI,.,NNI23:-(LLEIEJXALI.NNIJ-NNIXALI,NNI-*123)/ 
1።፣፤፡፤ ነ>!፤ 112752 
2=ህተ12 
NNI:-NILIJ-J; 
END; 
END; 
END: 
i =1+13 
END; 
END; 


PROCEDURE ONEATATIME; 


BEGIN 
IF NINITCTEJ=NICTIE THEN 
IE:=IE+1 
ELSE 
BEGIN 
NILIEJ:=NILIEJ+1; 
J:=0; 
REPEAT 


NN:zNIL[IEJ-J:; 
ALIE,NNJ:=LLIEIZALIE,NN-1I/ (NN-19; 
J:=J+1; 
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UNTIL NN=2; 


SUMMAT ION; 
SUBHI: 
END; 
END; 
PROCEDURE 511222 
BEGIN 
REL1:=0.0; 
IF TLITI=O THEN 
REL1:=1.0 
ESSE 
BEGIN 
FOR I:-1 TO K DO 
BEGIN 
SUM: 20.0; 
NKKs=NICII; 
FOR J:=1 TO NFK DO 
SOM: =SUMTALIT  JIIEXF((J-1)3ILN(TECIT1)):; 
ከ በስ ። 114 ር ረ፦(- 1 117,)1ቹጅቕ!!1131):5 
END; 
FORT STOD k DO 
الاجم‎ SRELI*PROLI J: 
END; 


1١ >21 1-1. 0 THEN 
RELETTIS-170 
EL SE 
PELLIT: SRELI: 
END; 


BIRMEEDURE PREPFORLOOFP: 
BEGIN 
FOR I:=1 TO MAX DO 
NINITCIIJ:=-1; 
FKK:=1; 
I:=1; 
fae 2.0 THEN 
BEGIN 
NINITL1]:=1; 
K:=1; 
END; 
WEILE LEI+123>0 DO 
BEGIN 
IA ALIAS THEN 
BEGIN 
NINITLKEKJ:=13 
K:=KKK; 
KKK:=KKK+1; 
s-141; 
«ን. ፡!፡ተ፤ ፡ር THEN 
NINITCKKKJ:=13; 
Ks =KKK 5 
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IF (X=-1) AND (Y=8) THEN PRINTVAR: =? 4? 
IF (X=-1) AND (Y=6) THEN PRINTVAR:="3" 
IF (X=-1) AND (Y=4) THEN PRINTVAR:="2" 
IF (X=-1) AND (Y=2) THEN PRINTVAR:=" 1" 

END; 


‘ma “40 We ‘48 


PROCEDURE GRAPHICS; 


BEGIN 
COUNT: =1; 
FOR Y:=20 DOWNTO 1 DO 
BEGIN 
FOR X:=-3 TO 76 DO 
BEGIN 
IF COUNT<=TMAX THEN 
BEGIN 
VAL : ZROUND (2OXREL E COUNT 3) ; 
IF (VAL=Y) AND (XZTLCOUNT3) THEN 
BEGIN 
PRINTVAR:=ASTERISK; 
COUNT : =COUNT +1; 
END 
EE SE 
PRINTVAR: =BLANK; 
END 
L'OE 
PRINTVAR: BLANK; 
YAXIS; 
WRITE (PRINTVAR) ; 
END; 
END; 
ሃ 5 =ርሂጋ2 
MORIA: =>=3 TO 7S DO 
BEGIN 


IF (X MOD 5 =0} THEN PRINTVAR:=LINE; 
IF COUNT<=TMAX THEN 
BEGIN 
VAL: =ROUND (20*RELC COUNT) ; 
IF (VAL=Y) AND (X-2TLCOUNT3) THEN 
BEGIN 
PRINTVAR: -ASTERISK; 
COUNT: =COUNT +15 
END 
ELSE 
TEA MOD 3 <> 90) THEN 
PRINTVAR: =DASH; 
END 
ELSE 
IF (MOD 5 <> 0) THEN 
PRINTVAR: =DASH; 
IF (X=-3) THEN PRINTVAR:=*0" ; 
IF (X=-2) THEN PRINTVAR:="."3 
IF (X=-1) THEN PRINTVAR:="0"3 





END 


ELSE 
BEGIN 
115 >1 
COUNT : z1; 
WHILE LEITI=LCII+13 DO 
REGIN 
COUNT: =COUNT+1: 
NINITLKEKKJ: COUNT ; 
BOR IE:-II TO 9 DO 
011-35: ح٠‎ 1 11+ 1[ 
LIMAX1:=-1: 
K:=KKK; 
END; 


KKK 3s =KKK +135 
Mo [Te 
IF (((COUNT>1) AND (LL£I+13<0) AND (LII3>0))) THEN 
BEGIN 
NINITLIKEEK 3:213 
K2=KKK; 
END; 
END; 
END; 
END; 


ERBBEDURE LOOP; 
BEGIN 
FOR IT:=1 TO TMAX DO 
BEGIN 
INITIALIZE; 
REPEAT 
ONEATATIME; 
UNTIL IE=K+1; 
= : 
NELETTIS=SRELLITIXPIBRI; 
ges OE | S2=RELL ITI; 
END: 
END; 


PROCEDURE 15: 

HEGIN 
IF (X=0) THEN PRINTVAR:=LINE; 
IF ((X=O) AND (Y=20)) AND (TL12=-0) THEN PRINTVAR: =ASTERISK; 
IF (X=-3) AND (Y MOD 2 = 0) THEN PRINTVAR:=? 0? ; 
IF (X=-3) AND (Y=20) THEN PRINTVAR: =’ 17; 
IF (X=-2) AND (Y MOD 2=0) THEN PRINTVAR:=”.”?; 
IF (X=-1) AND ((Y=20) OR (Y=0)) THEN PRINTVAR:="0"3; 
IF (X=-1) AND (Y=18) THEN PRINTVAR: =’ 9’ 
IF (X=-1) AND (Y=16) THEN PRINTVAR: =’ 8B’ 
IF (X=-1) AND (Y=14) THEN PRINTVAR: =’ 7’ 
IF (X=-1) AND (Y=12) THEN PRINTVAR: =’ 6’ 
IF (X=-1) AND (Y=10) THEN PRINTVAR: =’ 5’ 


De 49 28 we 8 
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WRITE(PRINTVAR)s 
END; 
WRITE (? O 5 10 13 20 29 30 23 40 45 
JO 55 60793 
WRITELN(' TIME <> JS 
END; 


PROCEDURE DECIDE; 
BEGIN 
WRITELN; WRITELN; 
WRITELN(’?DO YOU WANT TO CONTINUE WITH MIXED CONVOLUTIONS? 
IMARE "D" :"); 
WRITELN(’OR DO YOU WANT TO COMPUTE CONVOLUTIONS? 
INE "XP1" 2”): 
END; 


BEGIN 
INPUT; 
ENUTPUT1; 
GETDATAx 
BETPUT2: 
EEMPPIOR; 
COMPLAMB; 
FOR I:-1 TO MAX DO 
ሮ፣ 12=፦15 
FOR BR:=1 TO BRANCHES DO 
BEGIN 
s=13 
WHILE LAMBELSR,J1>=0 DO 
BEGIN 
IF i dB BR. 1120 TITHEN 
BEGIN 
ULJJ: =_5MBLER;/J1l1: 
J:=J+1: 
END 
ELSE 
BEGIN 
ELJISs-LAMBEBR,JE11; 
J:=J+1; 
END: 
END; 
SORT; 
REFE ORL OOP; 
EF: 
END; 
FOR I:=1 TO TMAX DO 
BEGIN 
SUMRELE IJ: =O; 
FOR BR:=1 TO BRANCHES DO 
SPRER 13s —=RELTOTEBRS TIFSUMREL EII: 
ሣው እ ኤእ ር ፡፡!|፡)ኀ፡11.!1|፣ርቄ27.,211112.,. ) 4 ”, 
SUMEEELII:S:6); 
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END; 

WRITELN; 
ATTENTION; 
WRITELN(*DO YOU ALSO WANT A GRAPHICAL PRESENTATION 
OF THE DATA?’ ): 
maine HIT "1" FOR YES OR "O" FOR NO : i 
READLN(CHGR) ; 
WRITELN; 
WRITELN; 
IF CHGR=1 THEN 

GRAPHICS; 
DECIDE: 

END. 
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